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Abstract 

Superintegrable systems in two- and three-dimensional spaces of constant curvature 
have been extensively studied. From these, superintegrable systems in conformally flat 
spaces can be constructed by Stackel transform. In this paper a method developed to es- 
tablish the superintegrability of the Tremblay-Turbiner-Winternitz system in two dimen- 
sions is extended to higher dimensions and a superintegrable system on a non-conformally- 
flat four-dimensional space is found. In doing so, curvature corrections to the correspond- 
ing classical potential are found to be necessary. It is found that some subalgebras of the 
symmetry algebra close polynomially. 

A maximally superintegrable quantum system on a n-dimensional manifold is an integrable 
Hamiltonian system of n mutually commuting differential operators and an additional n — 
1 differential operators so that the full In — 1 are algebraically independent and commute 
with one distinguished operator, the Hamiltonian, which we will take to have the form H = 
V 2 + V. In all previously known quantun superintegrable systems of this form with non- 
constant potential, V 2 is the the natural Laplacian of a constant curvature manifold. A Stackel 
transform can be used to construct systems on conformally-flat manifolds from systems of 
constant curvature manifolds [3j [8], but no systems have been previously exhibited on a non- 
conformally-flat manifold. 
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In [1] a classical superintegrable system on a non-conformally flat 4-dimensional space was found 
by generalising the two-dimensional Tremblay-Turbiner-Winternitz (TTW) system [T5] and here 
we will show that this system can be quantised in a way that preserves its superintegrability. 

The two-dimensional TTW system is integrable by virtue of a second order operator associated 
with its separability Its superintegrability has been demonstrated by constructing a symmetry 
from raising and lowering operators built out of special function recurrence relations that act 
on an eigenbasis of separated solutions |5J. A similar approach has been used to generate 
other families of superintegrable systems in two dimensions [5j [121 H31 EH] and in so doing, 
greatly expanded the list of superintegrable systems with higher order symmetries. Previous 
studies of higher order superintegrability have uncovered quantum superintegrable systems 
with no classical counterpart [JJ and the need to consider systems in higher dimensions with 
higher order symmetries has been highlighted recently by the use of higher order symmetries 
to determine the spectrum of a deformed Kepler-Coulomb system in three dimensions [T4"] . 

Here we extend the raising and lowering operator method used on the quantum TTW system 
to higher dimensions and construct sufficient additional algebraically independent operators 
to show that it is superintegrable. Furthermore, we find that some subalgebras of symmetry 
operators close polynomially as is common with previously known superintegrable systems. 

An interesting feature encountered below is that in order to construct the additional symmetries 
required for superintegrability the potential must be deformed by the addition of curvature 
terms that make the Hamiltonian conformally covariant. These terms are not simply the usual 
minimal choice for a conformally covariant Laplacian, namely, —1Z/6, where 1Z is the scalar 
curvature associated with the underlying metric, but also include an invariant constructed from 
the Weyl conformal curvature. 

The system considered below is 4-dimensional only so as to provide the simplest non-conformally- 
flat example. The procedure can be extended to higher dimensions with no greater difficulty 
and a number of other systems similar to the TTW system [5] could also be extended. 



1 The classical 4D non-conformally-flat system 

The Tremblay-Turbiner-Winternitz (TTW) system [T5] sparked great interest because it pro- 
vided an infinite family of superintegrable systems and examples of systems with arbitrarily 
high degree symmetries. The system in polar coordinates is given by 

U 2,2,1/2, A $2 

Httw =p r + ar + — I p g + ___ on _ m + 



cos 2 (k9) sin 2 (k9) 

The superintegrability of both classical [7] and quantum versions [2j were established for all 
positive rational values of the parameter k along with the polynomial closure of its symmetry 
algebra and the general approach was quickly extended to other families of systems in two 
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dimensions such as 



H = cosh 2 ip +p 2 v + 



a 8 7 

+ -4r— + 



cos 2 k(p sin kip sinh ip 



In the classical case, this approach as also been extended to higher dimensions and the 4- 
dimensional generalisation of the classical TTW system, 

H = L 1 = p 2 r + a r 2 + ^ (1) 
Ls = Pei + — I7T^T + 



cos 2 (Mi) sin^fc!^; 



r 2 , fa L 4 

3 ~ P92 + cos 2 (M 2 ) sin 2 (M 2 ) 

r 2 ^3 At 



cos 2 (M 3 ) sin 2 (M 3 )' 

was shown to be superintegrable for all positive rational k%, and &3 jl]. Here, the underlying 
manifold, on which this is a natural Hamiltonian system, has metric 

# = e ® e + ei <g> ei + e 2 ® e 2 + e 3 ® e 3 , (2) 

e = dr, ei = rd#i, e 2 = r sin(/c 1 ^ 1 )<i6' 2 , e 3 = rsin(/c 1 6'i) sin (/i^^) 6^3, 

and is no longer flat unless k% = k 2 = 1. Furthermore, in these coordinates, each component of 
the Weyl conformal tensor is a constant multiple of 
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r 2 sin (fci#i) 

and so the underlying manifold is only conformally flat when this quantity vanishes, that is, 
k\ = k 2 . In [4 J it was found that ([I]) is superintegrable with the 4 second order constants given 
above as well as an additional cubic and two quartic constants. 

A natural question to ask is whether there is a corresponding quantum system with the same 
or minimally modified potential. It is straightforward to check that 



H = d* + -d r -u 2 r 2 + ^ 



Li = $ + 2h cot{ki9 l )d 6l + — f[ r + — 4jT a , 
1 cos^fc^i) sin (kiOi) 

L 2 = dg 2 + k 2 cot(k 2 6 2 )de 2 



cos 2 (k 2 6 2 ) sin (k 2 d 



2j 



T «2 P3 P4 

^3 = dg 3 + „„ - x + 



cos 2 {k 3 9 3 ) sin 2 (k 3 6 3 
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are four mutually commuting differential operators and H is a Hamiltonian of the form H 
V 2 + Vo where V 2 is the Laplacian on a 4-dimensional manifold with metric and 

T/ 2 M 01 M , 03 



r 2 cos 2 (/ci^i) r 2 sin 2 (/ci6>i) cos 2 (/c 2 # 2 ) ?" 2 sin 2 (fci0i) sin 2 ^^) cos 2 (£;3#3) 

+ & ( 3 ) 

r 2 sin 2 (A; 1 ^i) sin 2 (M 2 ) trin 2 (Jfeg0g) ' V ; 

However, it is not a simple matter to quantise the additional classical constants found in jl], 
and so to investigate whether this system remains superintegrable in the quantum case we 
attempt to adapt the raising and lowering operator methods from [5]. 

In the two-dimensional examples, with parameter k, the essence of the method is that solutions 
can be found by separation of variables with the separated eigenfunctions enumerated by posi- 
tive integers no and n\. The energy eigenvalue of each separated solution depended only on the 
combination no + kn\. Differential operators were then constructed to raise or lower no or n\ 
by integer amounts and so that for rational k, compositions of these operators could be found 
that left n + k\n\ unchanged and hence preserved the energy. While these additional operators 
were constructed to act on a separated eigenbasis, they were found to in fact be expressible as 
differential operators. It should be noted here that the linear dependence of the energy on the 
quantum numbers no and n\ is crucial to the method and maintaining this below leads to the 
need for quantum corrections to the potential. 

In order to extended this approach to our present example we must first solve the system by 
separation of variables and so we postulate a solution to = E^fr of the form 

with 

L 3 ^3 = £ 3 ^3, L 2 * 2 ^3 = 4*2*3, il*l*2*3 = 4*1*2*3- 

While the solution of the separated equations is unremarkable, the details are written out at 
length so as to expose the point at which the quantum corrections (jHJ) and ffT2l) to the potential 
become necessary. 



U '\V) + ( ^-2— + ^77 + (2« + a + + l) 2 ) u(y) = 0, 



We find that each angular equation is, up to a gauge scaling, of the form 

1 « 2 

snT y cos 2 y 

which has solution 

u{y) = (sin y) a+ ^ (cos yf + ^P^ (cos 2y) 
where Pi Q,/3 Vx) is a Jacobi function HU1. 



Starting with the 9 3 equation, we make the replacements 

03 = k 2 3 ( J -a 2 J , /3 4 = fc 2 Q-a 3 
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and find the separated equation is 

y cos 2 (k 3 9 3 ) sm 2 (k 3 6 3 ) J 

which has solutions 

KC(0s) = (sin(A;3^)) a3+ ^cos(A;3^)) a4+ ^£ 3 ' a4) (cos(2A;3^)) 

with eigenvalues 

L^tC{e 3 ) = £^ a 3 %T(0 3 ), i 3 = -k 2 3 (2n 3 + a 3 + a, + l) 2 . 
The separated equation L 2 ^ 2 (9 2 ) = ^2^2(^2) is now 

9W 2 ) + k 2 cot(k 2 9 2 )%(9 2 ) + ( f 2 + , * 3 ) ^2(82) = £2^2(62) 

\cos z {k 2 2 ) $vcf(k 2 9 2 ) J 

which we transform with 

V 2 (9 2 ) = (sm(k 2 9 2 ))-h 2 (6 2 ) 
to absorb the first derivative term to give 



+ ( — ^— + $+j£L + 1 -ki-i 2 ) l mm = 

Vcos^A^J sin (« 2 02j 4 
and we make the replacements 



1 l\-<4)> £ 3 + \ki = kir--Ai h 



which when combined with (j2J) gives 



k 3 

A 2 = -r{2n 3 + a 3 + a 4 + 1). 
k 2 



The separated 6 2 equation becomes 



where 



cos z {k 2 9 2 ) sm (k 2 9 2 ) 4 



k 2 

-j- - £ 2 = k 2 2 (2n 2 + a 2 + A 2 + l) 2 



and has solution 

*2nf(^) = (sin(A;2^)) A2 (cos(A; 2 ^)) a2+ 5P^^)(cos(2A; 2 2 )). 



The separated 61 equation L 1 \& 1 (6' 1 ) = £ 1 \& 1 (6' 1 ) is 

+ 2h oot(Mi)*i(fli) + f Ynr~a~\ + . Jhn ) = WO, ( 7 ) 

\cos z (ki^i) sm (fci 6*1)/ 

which we transform with 

= (sin(Mi)rV(#i) 

to absorb the first order term to give 

VWO + f iTbl + • , + *? - V*(*i) = (8) 

\co^{kiOi) sin (kiOi) J 

and we make the replacements 
Combining this with (IB]) gives 



A = W^(l-|)+|(2- 2 + a 2 + A 2 + l) 2 . 

This does not have the same form as (jSJ) and will not lead to an energy eigenvalue that depends 
linearly on n 2 . Hence, we instead propose an additional quantum correction in the potential of 

r z sin ykiUi) 

which in turn leads to a modified (jSJ), 

+ (^^r + i2 + fi:f 2) + *? - 4) v»i(flx) = 0. 

\cos^(rei0i) sm (fcifc'i) / 
Now, making the replacements 

ft = *? a?) , ^ + \{kl-kl) = kl{^-A\ 

gives 

A 1 = ^(2n 2 + A 2 + a 2 + l). (10) 
The separated 6\ equation is now 



where 

k\ - t x = kl(2m + A 1 + ax + l) 2 

and has solutions 

*&f(0i) = (sin(A; 1 ^ 1 )) A ^5( cos (A; 1 ^ 1 ))^+5p ? (^^)( cos (2A; 1 1 )). 
Finally, the separated radial equation is 

i?*o(r) = 9 2 ^ (r) + ^9 r *o(0 + f-wV + ^ * ( r ) = E^ (r). (11) 

In a similar way to above, in order that E depend linearly on ni, we propose the addition of a 
quantum correction to the potential of 

V,= l -± (12) 

which leads to a modified version of (jlip . 

H%(r) = d 2 r Mr) + ^Mr) + f^V + £l ~ ^ — \ Mr) = ^o(r). 
We remove the first order terms with the transformation 

3 

^o(f) = r~z?po(r) 



to give 



Now, 



has solution 



%Mr) + ( -"V + I ^ + - E ) Mr) = 0. 



- - A 2 

u"(x) + f -x 2 + 4 - + 4n + 2A + 2 ) w(x) = 



where Ln (x) is a Laguerre function |10j . 

We needed A\ = k\ — i\ and we already have k\— 1\ — k 2 (2rii + A x + a\ + I) 2 so 

*fio(r) = ^-^"^(wr 2 ), (13) 

where the mulitplicative factor of w Ao/2 is chosen for later convenience, and 

A = fci(2m + ai + Ai + l), E = -w(4n + 2A + 2). (14) 

Note that with the quantum deformation ( TT2|) the relationship of to n\ is similar to that 
seen in (J5J) and ffTOj) . 



Now, putting together (jSJ), f lTOj) and f JT4|) we find 

= -2w(2n + 2/^ni + 2k 2 n 2 + 2A; 3 n 3 + fciai + k 2 a 2 + ^30.3 + h a 4 + h + k 2 + k 3 + 1) (15) 
for a solution of the form 

^no,ni,n 2 ,n3 *0,rio v J * l,m \ u l)^2,n 2 \ w 2J*3,n 3 \ W 3J- 

With the quantum corrections Qj and ffT2l added to the potential we have the following set of 
mutually commuting differential operators. 



Li = ^ 1 + 2fcicot(Mi)9 fll + ij^—rr + . 2 , , . : + 



cos 2 (Mi) sin 2 (Mi) 4sin 2 (M 



^2 = <\, + &2 COt(fc 2 2 )^ 2 H 2 / ? fl x + 



cos 2 (k 2 9 2 ) sin 2 (k 2 6 2 ) 



T - r) 2 1 ^ 1 & 



cos 2 (k 3 6 3 ) sm 2 (k 3 3 )' 

= E^fr remains separable with the additional terms. In the following, we use these redefined 
H and Li. 

For the metric (jSJ), the scalar curvature is 

K - 6 k 2 ( 6 2 ^ 2A;2 

r 2 1 \r 2 r 2 sin 2 (ki9 1 ) ) r 2 sin 2 (£;i6>i) ' 

and with Weyl conformal tensor W a b c d, if we define 

w = Vzw abcd w abcd = r^3- v 

rsm {kiVi) 

then 

# = V 3 + Vb + Vi + V r 2 = v 2 + y -^-^->v. 

o 24 

Note that V 2 + V\ + V2 is a conformally covariant Laplacian and the metric g is conformally 
flat if and only if k\ = k 2 . 



2 Raising and lowering operators 

Our aim is now to use special function identities to raise and lower the n, while preserving E 
and produce new operators commuting with H. 
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Using differential identities for Laguerre functions [TO] we construct the operators that act on 
the radial part of the separated solutions, 

K+ n A ; = ^^ r + (27^ + Ao + lV + i^^ 



Ko n A ° = 1 -^d r + (2n + A + l)u+ 1 —^. 



These raise or lower no by 1 while simultaneously lowering or raising and Aq by 2, that is, 

Ktn A °<l = -2a;(n + l)(no + Ao)< +i: 

Note that the constant multiplicative factor of tu A °^ 2 in ffT3|) was chosen so that both of these 
have a factor of u on the right hand side. 

For the angular functions, we can use Jacobi function identities [TU] to make operators that 
raise and lower n alone, 

J: = - {N + 1) 2 f ,(m 9, - \ ((N + 1) (N + 1 - c - i) cos(2^) 
- (A^+ l)(c-d) + a 2 -6 2 ), 

Jn = <JV " l ^ ke h e -\((N- 1)(N -l + c + d ) cos(2^) 
+ ( N - l)( c _ d) + a 2 - b 2 Y 

where N = 2n + a + b+ l and their action on 

Q(a,b) = sin a+c(k0) COS b+d (k6)Pl a > b) COs{2k6) 

is given by 

= -2(n + l)(n + a + Hl)eS, 



J-Q^ = -2(n + a)(n + b)e 



(a,b) 
n-1 ' 



The operators above are essentially those used in [5] and the analysis used to show superin- 
tegrability for the TTW system immediately carries over the the current example and so we 
obtain a symmetry operator by raising and lowering the functions associated with r and Q\. 

Notice that the operators above that raise or lower the radial eigenfunctions change both uq and 
Aq. In order to extend this approach and construct a symmetry operator by raising and lowering 
the 9\ and 8 2 functions we will need operators with a similar effect on these functions. This 
can be achieved using Jacobi function identities that raise and lower n and a simultaneously 
when applied to 

Q { n a ' b) = sm a+c (k6) cos b+d {kd)P^ b) cos(2k9). 
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We find 



with action, 



Ky = - (1 k 0) d e - 2(n(n + a + b + 1) + a(a + 6)) 

-(l-a)(fl + c + 6 + d)- (1 - a)(a " c) 



sin 2 (^) 

^ „ (1 + a) cos(k9) n n , 

- -4iWr^- 2n(n+a+6+1) 

-(l + a)(a + c + 6 + rf)+ (1 + a)(a + c) 



sin 2 (fc0) 

^ n+ a M) = 2 (n + l)(n + a)efc 2 '*', 
iT- a ei a ' b) = 2(n + a + 6+l)(n + 6)0l a+ 1 2 ' b) . 



3 Constructing the symmetries 

For fci = pi/qi with gcd(pi, gi) = 1 the operator 

_ A)+2(pi-l) K - A T+ T+ 

"1 ~~ xv no-(pi-l) " ' -"-0 n J l ni+gi-l ' " J l m 



Pi terms qi terms 

has the effect on a basis function of 

n ->n -pi, ™i ~> ™i + A ->Ao + 2pi, 

and so 

£ = -2w(2n + 2fcini H ) ->• -2w(2(n - pi) + 2fci(ni + gi) + 

= -2w(2n + 2£;ini H ), 

that is, E is unchanged. A similar lowering operator is 

W~ — K + A o-2(pi-l) K + A j- j- 
"1 ^0 no+(pi-l) n ^1 ni-(gi-l) 17 1 m ' 

N v ' V v ' 

Pi terms gi terms 

which also leaves E unchanged, has the effect on a basis function of 

no -> ™o +Pi, ni -> - ?i, A A - 2p x . 
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Explicitly, the action of on a basis function is 

Sf*&,*&x Bl = (-2) pl ^H^i + l) 9l K + A + ai + l) 9l <^^ gi ' 

S^nXnf = (-2)*V*(-m - A 1 ) qi (-n 1 - ai ) qi (n + Pl ) Pl (n + A,) P X%~l P P \^tT- qi - 

This is exactly like the TTW raising and lowering operators from [5]. 

For fc 2 /fci = P2IQ.2 with gcd(p 2 , Q2) = 1 the operator 

^+ _ K - a 1+ 2( P2 -i) K - Al j+ ...7+ riRi 

"2 — ^1 ni _(p a _l) ^1 ni J 2 n 2 +q 2 ~l J 2 n 2 
v v /V v ' 

p 2 terms <?2 terms 

has the effect on a basis function of 

rix ->■ ni — £>2, n 2 n 2 + g 2 , Ai -»■ Ai + 2p 2 

and so 

E = -2co(2n + 2fcini + 2/c 2 n 2 + •••)-► -2w(2n + 2fci(na - p 2 ) + 2A; 2 (n 2 + g 2 ) + ■ ■ • ) 

= -2w(2n + 2fci?ii + 2/c 2 n 2 + •■■), 
that is, .E is unchanged. A similar lowering operator is 

"2 - ^1 „ 1+ ( p2 _l) -fM m J 2 „ 2 _( g2 _l) J 2 n 2 • l ly J 
v v ' v v ' 

P2 terms g 2 terms 

This is similar to the two-dimensional TTW procedure, but different operators are required. 

The operators given so far are only well defined on the separated basis functions and they 
contain the quantum numbers in their defintions. We now must show that we can construct 
pure differential operators. The argument is only sketched here as the details are essentially 
the same as those in [5]. 

The transformation n\ — > —n\ — A\ — a\ — 1 while holding E constant has the effect of changing 
the sign of A . It is then straightforward to check from the explicit expressions for the operators 
that + 

Lt = Et + E7 and L7 = k ^ 1 - " 1 

A 

2 A 2 



are polynomials in A$ and A\. Since 



A\ = k\ - h and A\ = ^ 



we can replace E, A^ and A\ with second order differential operators where ever they appear 
in these expressions. 
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Similarly, the transformation n 2 — >■ —n 2 — A 2 — a 2 — 1 while holding L\ constant has the effect 
of changing the sign of A\. It is then straightforward to check from the explicit expressions for 
the operators that 

L+=~+ + ~ 2 and L 2 =k 2 ^ 2 (20) 
are polynomials in L 1; A\ and A\. Since 

A\ = ^~r 2 and A\ = --^ 

we can replace L\, A\ and A| with a second order differential operators where ever they appear 
in these expressions. 

In a similar way, we can define and Lf with the label replacements 1 — > 2 and 2 — )• 3 in 
(fT8]h ( TT9|) and ( 120]) and show that they are in also differential operators. 

It is clear from the construction that {H, L\, Lf, L 2 , Lf, L3, Lf} forms an algebraically inde- 
penent set of differential operators and hence the system is superintegrable. 



4 The symmetry algebra 

A common feature of superintegrable systems is a polynomially closed symmetry algebra. By 
direct calcuation, we find some polynomially closed subalgebras of the symmetry algebra. 

Adapting the an argument from [5] we find that, for % = 1, 2, 3, 

P^(L t _ 1 ,L l ,A 2 ) = E~Ef + EfE- and ^(Lw, L h A 2 ) = h~ + ~ 



are differential operators that are polynomial in their arguments. 

By comparing the action of brackets of the operators {Li, Lf, L~} with symmetrised products 
of the operators we find the following explicit identities for i — 1, 2, 3. 

[Li,L7] = -Ak 2 q 2 L~ - Aa^qiLf 

[L h Lf] = 2q i {L i ,Lr}-4k 2 q i Lf + 4k 2 q^Lr + 8qfk 2 Lr 
[Lf,L~] = 2q l {Lrf-2Pt\L l ^L l ,Al) 

and 

{L t , L7, L'} + 2k 2 (Uq 2 - 3a,)(Lr) 2 + 6k 2 (Lf) 2 + 6Afo{L+ L;} - UkfP^ + AkfaP^ = 0, 



where a% = 1, a 2 = 1/4 and 03 = 0. These hold as operator identities on general functions. 

'4,Lf\ = [Lf 



Furthermore, [Lj, Lf] = for i ^ j and [Lf, Lf] = [Lf, Lf] = for \i - j\ > I. 
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As was found for the TTW operators, the symmetries constructed from raising and lowering 
operators are not necessarily of minimal order The same technique for finding lower order 
operators can be used for the current system. 

For example, starting from Lf we look for satisfying 

[L 1 ,Mf]=Lf, 
[H, Mf] = 0, [L 2 , Mf] = 0, [L 3 , Mf) = 0. 

Find Mf is 

M~- 1 ( Lr | — ^ | S ^ L ^ (21) 

where Si(H, L 2 ) is a polynomial in H and L 2 that can be determined by the methods used in 
0- 



5 An example 

Explicit computations can be performed for particular choices of the fcj. For example, with 
ki, k 2 , k 3 = 2, 1, 1, the operator Lf is 6 th and L\ is the 5 th order operator, 

/sin(4#0 „ 3 cos(40i) + 2 - a? 5 s in(40O _ (6 cos(40i) + 5 - 4a?) \ _ 

~ y " °r°0i + " 6 r 4^ OOl J 

sin(40O 2 13sin(4^)_ (4 cos(4^) + 1) 2 13 + 27cos(4^) - 4a? 

+ 73 C'rC'O! =-5 »r + °r 



4r 2 r -1 4r 3 ' -1 2r 2 r 2r 3 

5sin(40i) 25cos(45i) + 20- 12a? \ 2 
~ *i 2^ ) A ° 

llsin(4^) (ll-8a? + 14cos(4^)) 23sin(4^) 

+ ^2 + ^ °V 4^ °r°0i 

26 0)8(450 + 23- 20a? _ 21sin(40i)_ 3(lOcos(40i) +7 - 4a?) 

O r — ; Of)-, — 



2r 3 4r 4 1 2r 4 

with the replacements E — > H, Aq — > k\ — L\ and A\ — > {k\ — 4L 2 )/(4/c?). 
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A 4 th order operator can also be constructed using (|2ip and in this case, 

qmr , (H 2 — Au)(A\ — q|) 
b 1 (H,L 2 ) = — 

and Mf is a polynomial in E and even powers of A and Ai. 



6 Conclusion 



The methods developed in [5] have been extended to demonstrate the superintegrability of 
a 4-dimensional quantum Hamiltonian system on a non-conformally-flat space. The system 
discussed is a quantisation of a previously described classical superintegrable system and in 
order to maintain superintegrability in the quantisation, correction terms were required to be 
added to the potential. These correction terms make the Hamiltonian conformally covariant, 
but are not the usual minimal conformally covariant correction of —71/6 as they depend on 
the conformal curvature. The 3-dimensional analogue of this system also requires the addition 
of the term —71/8 to maintain superintegrability and this too gives a conformally covariant 
Hamiltonian. 

While many previously known superintegrable systems possess a polynomially closed symmetry 
algebra, here we have only found some polynomially closed subalgebras. An investigation 
of a closely related classical system found that in general the symmetry algebra will close 
rationally rather than polynomially [6]. It seems reasonable to conjecture that, except in some 
special cases, the symmetry algebra of the 4-dimensional system considered here does not close 
polynomially, but rather it obeys an appropriate quantum analogue of rational closure. 
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